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O ' 1 Abstract 

(N 

Formation of luminous matter in the Universe is a complicated process, which includes many processes and 
components. It is the vastly different scales involved in the process (from star formation on few parsec 
scales to galaxy clusters and superclusters on megaparsecs scales) and numerous ill-understood processes, 
QQ \ which make the whole field a maze of unsolved, but exciting problems. We present new approximations 

CN ' for numerical treatment of multiphase ISM forming stars. The approximations were tested and calibrated 

CN \ using iV-body+fluid numerical simulations. We specifically target issues related with effects of unresolved 

l^^_^ ' lumpinesses of the gas. We show that the degree of freedom is much smaller than naively expected because of 

^j. \ self-regulating nature of the process of global star formation. One of the problems of numerical simulations 

is related with the uncertainties of approximation of the supernovae (SN) feedback. It is often assumed that 
the feedback is mainly due to momentum transferred by SN in to the ISM. We argue that this may not be 
true. We present a realistic example of gas actively forming stars with short cooling time for which the SF 
O \ feedback is important, but the kinetic energy of the gas motion due to SN is negligibly small as compared 

-4— » ' with the thermal energy of gas. 
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2 Introduction 



Numerical A^-body+hydro simulations are important tools for making theoretical predictions for galaxy 
C^ ' formation in the expanding Universe. Unfortunately, the simulations are mired by two basic problems, i) 

Results on scales of interest (kpc and larger) crucially depend on much smaller (pc) scales. This is very 
different from A'^-body problem where the situation on large scales does not depend on what happens on 
small scales. The reason why small scales are important lies in the nature of the cooling and heating of the 
gas. Small-scale lumpiness changes the gas cooling rate. The gas heating is defined by how the star formation 
proceeds on very small scales, not simply by the average star formation rate, ii) There is no theory of star 
formation (SF). Observations regarding SF are very useful, but they do not cover all important situations. 
Nevertheless, the situation is not that grim as it looks. Not all details of the small-scale processes are 
really important. Recently, Yepes et al.(1997)(l] presented realistic simulations of the galaxy formation in 
which different parameters (e.g., cooling rates, resolution, SN feedback) were changed dramatically with the 
goal to find out how the final results are sensitive on the parameters. The final global results (for example, 
total luminosities and SFR) were remarkably stable when the cooling rates changed by a factor of ten and 
when the resolution changed by a factor of two. The results for small galaxies were very sensitive to the 



feedback parameter, while large galaxies were insensitive even to the SN feedback. Some of the dependencies 
are easy to understand, but some are more difficult. For example, if the cooling time is very short (as 
typically is the case), the cooling time is not important. Another reason for stability lies in the nonlinear 
response of the ISM to the SF. The ISM is typically found in a self-regulating regime when the rate of infall 
of a fresh gas from outside is defined by the heating of the gas by newly formated gas. If the infall is large, 
the SF increases (because the gas cools and becomes available for SF), this results in larger heating of the gas 
by young stars and SN, which stops the infall. This self-regulating evolution of the ISM results in effective 
"canceling" of the number of free parameters ~ results are sensitive only to a small number of combinations 
of the parameters. 

The following example illustrates the idea. Let x and y be two variables, which describe the situation 
(e.g., density and temperature of gas). They depend one on the other in a simple, but nonlinear way, which 
involves four parameters oi, 02, 61, ^2- The evolution of the system is defined by two differential equations. 
The system has the same basic structure as a real system: ^ = oix — 612/^, -^ = 02?/ — b2X^ . The system 
depends on four parameters 0^,6^ and on two initial conditions: xq, 2/0 ■ Thus, six "free parameters" define 
the evolution of the system in general. With six free parameters one naively would expect that any final 
configuration is possible. But in the regime of self-regulating evolution, ^ ~ ^ ~ 0, and the solution 
depends only on two parameters a = (ai/61)^/^ and /3 = (02/^2)^^'^, and it does not depend on initial values 
of X and y. 

3 Cooling of gas with multiple phases 

The equations describing the evolution of a two-phase medium of hot gas, which emits radiation, and cold 
clouds (Tc — 10"*-?^), which effectively have stopped cooling, but are capable of producing luminous matter 
("stars"), can be written in the form pi: 



den P , ., ,, 2^(^/0 insi 
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Ph^^^[esN-A{e,,^e,)]^apl^^, esN = Werg/22MQ, (1) 



^--^^P^.4^ + ^, /3-0.12, (2) 

at Eh — tc n^'mjj r* 

where indices h, c, and gas refer to the hot, cold, and the total gas components; kT = (7 — l)mH^im^■ Free 
parameters a and C were calibrated using numerical simulations: a — 0.95, C = 2. The feedback parameter 
must be in the range A — bQ — 200. In order to test our approximations for cooling of multiphase gas 
we run few numerical simulations of evolution of initially slightly inhomogcncous gas without gravity and 
without SF. Simulation box is small: 500pc-3kpc. For a real cosmological simulation the whole box would 
be one resolution element. We use 64"^ PPM code to run the simulations. Initial conditions were either small 
random gaussian fluctuations with RMS=0.15 or a simple sine-wave perturbation 5p oc sin(a;) sin(j/) sin(z). 
Two combinations of the size of the box and the gas density were chosen; i) gas cools very fast (no motion 
of gas) or ii) the gas cools slowly, and it has time to produce motion of gas across the box. Results are 
presented in Figure 1. 

We also used the initial conditions to make runs with the SF and feedback included. In all cases we 
observed that after some initial evolution the system settles in a regime of steady "burning" of gas into 
luminous matter ( "stars" ) . In that regime the energy produced by stars is almost equal to the energy 



Figure 1: Cooling of gas in the regime of thermal instability. Full curves are for gas with initial random 
distribution of density fluctuations with the rms a = 0.15. Initially the gas is in the pressure balance and 
T = lO^K. The dot-dashed curves are for gas with one sine-wave perturbation of density with the same 
rms fluctuation as for the full curves. As the gas cools the thermal instability develops. It results in the 
formation of very dense cold lumps and hotter gas with very low density. The long-dashed curves show 
evolution of a homogeneous gas with the same initial average density and temperature. The dotted curves 
show results for our multiphase treatment. It clearly provides much better approximation as compared with 
the homogeneous gas. 
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radiated by the gas with the kinetic energy of the gas always being much smaller than the thermal gas 
energy. In the regime of the self- regulating SF the gas temperature appears to be a constant, which can be 
found analytically: 

_ (C - a)PesN , (C - a) esN ^sn ..^ 

^'^ - C(l - /3 - A/3) + '" ^ C A ^ A ^ ' 

Figure 2 presents an example of the evolution of the multiphase gas with the supernovae feedback. The 
SF was delayed by one cooling time in order to increase the fraction of energy in the form of gas motion. 
Even in this extreme case, the kinetic energy was much smaller than the thermal energy of the gas. After 
2 X 10 yrs the system was already in the regime of self-regulating SF. The total energy released by forming 
stars was almost equal to the energy radiated by the gas with the kinetic energy being much smaller than 
each of the energies. The temperature of the gas was close to the value predicted by equation (4). Note that 
vastly different initial conditions (one sine-wave or 1/4 million independent fluctuations) resulted in very 
similar final states of the gas and almost identical rates of conversion of gas to stars. 
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Figure 2: Evolution of average gas properties in the case of active star formation and the supernovae feedback. 
3D PPM hydro code with 50pc resolution was used. The star formation was allowed only after lO^yrs. The 
effective time-scale for SF t*e// = 1-5 x lO^yrs was much longer than the minimum t^, — 5 x lO^yrs imposed 
by the code. 
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full = initial gaussian fluctuations, ct = 0.15 
dash = one wave, initial amplt=0.30 
Box = 3.2kpc Ax=50pc 
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